X-643-65-430 


gpo PRICE $ 
CFSTI PRICE(S) $ 

Hard copy (HC) 
Microfiche (MF) 

ff 653 July 65 


NASA 


1 


£ r p r- r <i r.J 

hj A" 


553*10 


INJECTION CONDITIONS 
FOR 

LUNAR TRAJECTORIES 


BY 

RONALD KOLENKIEWICZ 
WILLIAM PUTNEY 

N66-17262 



(ACCESSION NUMBER) 

(THRU) 

/ 

(PAGES) 

(CODE) 

3 Q 

(NASA CR OR TMX OR AD NUMBER) 

NOVEMBER 1965 

("CATEGORY) 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 


INJECTION CONDITIONS FOR LUNAR TRAJECTORIES 

by 

Ronald Kolenkiewicz 
William Putney 


November 1965 


Goddard Space Flight Center 
Greenbelt, Maryland 





INJECTION CONDITIONS FOR LUNAR TRAJECTORIES 

ABSTRACT 

This paper reviews the mechanics of Earth to Moon trajectories as 
affected by geometric considerations and booster capabilities. It formu- 
lates and describes the equations for a computing procedure which, using 
the two body equations of motion, provides approximate initial injection 
conditions near the earth for either a fixed time of arrival at the Moon or 
a fixed time of flight to the Moon. Consideration of multiple circular 
parking orbits, arbitrary injection elevation angles, arbitrary launch site, 
and booster burning characteristics are also taken into account. The 
trajectory followed is an elliptical Earth to Moon transfer trajectory, with 
respect to the Earth, which intersects the Moon before apo-apsis. A 
digital computer program having the equations programmed in FORTRAN 
II is available. 
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INJECTION CONDITIONS FOR LUNAR TRAJECTORIES 


I. INTRODUCTION 

A calculating procedure, using two body equations of motion, has been de- 
veloped which provides approximate injection conditions for earth to moon traj- 
ectories. The resulting injection conditions are intended for use as a first guess 
in generating precision trajectories. When the results obtained are used in a 
precision digital computer program (for example ITEM, see Reference 1) they 
typically produce a lunar impact. They can be further adjusted by iterating in 
the precision program to provide a specified miss distance for such missions 
as a lunar orbiter. 

The procedure presented builds on and extends the work of References 2 to 
4. In Reference 2 a study was made to determine the effects on lunar trajectories 
of some typical geometric constraints. This study concludes that application of 
the various constraints seriously restricts the allowable launch times during the 
month and day for direct-ascent launch; whereas, less serious restrictions re- 
sult for the parking orbit launches. This reference also presents some of the 
equations for the calculation of geometrical parameters involved. Reference 3 
further justifies the parking orbit type of trajectory. It also presents some of 
the equations necessary for matching the powered phases of the trajectory to the 
geometrical constraints. Reference 4 gives some general discussion on lunar 
trajectories. 

The present paper adds the consideration of multiple circular parking orbits, 
arbitrary injection elevation (or flight path) angle, and booster burning charac- 
teristics. Injection conditions can be determined for either fixed time of arrival 
at the moon or a fixed flight time to the moon. The mechanics of earth to moon 
trajectories as affected by geometrical considerations and booster capabilities 
are reviewed, and the working equations which are used in a computer program 
are formulated and described. 
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II. NOMENCLATURE 


x, y, z 


I. J> K 


D, 


r 

m 


my 


t 

m 


D 


m 


i 

A 

a 

S 

w 


w 


w 


w 

z 


V 

¥ 


L 


Components of geocentric equatorial coordinate system 
(Figure 2). 

Unit vectors along geocentric equatorial coordinate system 
(Figure 2). 

Unit vector to launch site. 

Launch time, hours. 

Number of days from January 0.0 U.T. 1960 to 0.0 U.T. on 
the day of launch. 

Unit vector to the Moon. 

Components of the unit position vector to the Moon. 

Time of impact at the Moon, hours. 

Number of days from January 0.0 U.T. 1960 to 0.0 U.T. on 
the day of arrival at Moon. 

Inclination. 

Azimuth (geocentric). 

Right ascension. 

Declination. 

Unit vector normal to the orbit plane. 

Components of the unit vector normal to the orbit plane 
(Figure 2). 

Angle between the K and wx r l vectors, see Figure 2. 

Mean longitude of Sun's apparent motion around the Earth 
defined by Equation (14), radians. 

Longitude (terrestrial). 
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Absolute rate of spin of the Earth. 


Total time of flight, defined by Equation (33), hours. 

Distance from the center of the Earth to the Moon. 

First boost burning time, ( t t - t ; ), hours. 

First boost angle, Earth centered angle between R ; and R t . 
Time of injection into parking orbit, hours. 

Time in parking orbit, (t 2 - t t ), hours. 

Reciprocal of the parking orbit rate, defined in Equation (18). 
Parking orbit angle, Earth centered angle between R t and R 2 . 
Earth’s gravitational parameter. 

Distance from the center of the Earth to the parking orbit 
injection. 

Angular distance (Figure 5). 

Second boost angle, Earth centered angle between R 2 and R i . 

Translunar trajectory angle, Earth centered angle between R ; 
and R . 

m 

Second boost burning time, ( t. - t 2 ), hours. 

Time of injection into translunar trajectory, hours. 

Time at beginning of second boost, hours. 

Time in translunar trajectory, ( t m - t. ), hours. 

Semimajor axis. 

Mean anomaly. 


Eccentricity. 


<<? 


E Eccentric anomaly. 

p Semilatus rectum of an ellipse, a (1 - e 2 ) . 

v True anomaly. 

R Distance from the center of the Earth to injection into trans- 

i 

lunar trajectory. 

R Ratio of injection distance, Ri, to moon's distance, R m . 

y Injection elevation angle (See Figure 5). 

V. Injection velocity. 

V p Parabolic velocity. 

Ratio of injection velocity, V j, to parabolic velocity, V p . 
R^g) Rotation matrix, defined by Equation (42). 

R 2 (8) Rotation matrix, defined by Equation (43). 

R 3 (0) Rotation matrix, defined by Equation (44). 

Y Year of lunar impact. 

m 

R ; Distance from the center of the Earth to launch site. 

h po Parking orbit altitude, 

h. Injection altitude. 

V c Circular velocity. 

V m Lunar impact velocity. 

R g Earth's radius (taken as 6378.165 km). 

V cm Circular velocity at lunar distance. 

T tA Equivalent to T t defined by Equation (15). 
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T tB Equivalent to T t defined by Equation (33). 

AT t Total time difference, defined by Equation (38). 

Vj Minimum value of the velocity ratio, defined by Equation (34). 

V 2 Maximum velocity ratio, defined by Equation (35). 

(AT^j Total time difference, Equation (38) , for V , see Figure 6. 

(AT t ),, Total time difference, Equation (38) , for V 2 , see Figure 6. 

V k 1 Second and succeeding assumed velocity ratios, defined by 

Equation (39). 

Final value of assumed velocity ratio, see Figure 6. 

S Parameter defined by Equation (37). 

R mx x component of the moon position vector. 

R my y component of the moon position vector. 

R mz z component of the moon position vector. 

SUBSCRIPTS 

/ Launch. 

1 Start of parking orbit. 

2 End of parking orbit, 

i Injection. 

m Moon, 

n The orbit number. 
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III. GENERAL DISCUSSION OF THE PROBLEM 


A. Requirements for Translunar Trajectory 

In general, there are two requirements that must be met for launching a 
vehicle from the surface of the Earth to the vicinity of the Moon with a minimal 
amount of energy. 

The first is that the vehicle be launched into a plane that is common to the 
launch site at launch time and the Moon at impact time. To launch in a plane 
other than this necessitates changing the plane of the trajectory after launch. 

This entails a reduction in the payload that will arrive in the vicinity of the Moon 
since additional fuel is required to accomplish this maneuver. 

The second requirement is that a specific angular relationship must exist 
in the plane between the center of Earth and launch site and the line of centers 
of Earth and Moon. This comes about from the fact that there are a number of 
restrictions on various lunar missions such as booster capabilities, allowable 
time in parking orbit, and the time requirements on the mission. 

The following sections will discuss the type of trajectory that must be fol- 
lowed to satisfy requirements on both launch plane and angular travel in the plane. 

B. Launch Considerations 

In general, for any given launch azimuth there are two times a day a vehicle 
may be launched from the Earth into a trajectory that would take it to the Moon. 
The exceptions occur only when the absolute value of the Moon's declination is 
greater or equal to the inclination of the trajectory plane; then, either none or 
at most one launch per day will be possible. 

C. Boost Considerations 


There are two possible methods for leaving the surface of the Earth and 
going to the vicinity of the Moon. The first is lunar injection from a direct 
launch; the second is lunar injection from a parking orbit. These trajectories 
are shown in Figure 1. Both the lunar injection from a direct launch as well as 
from a parking orbit have been discussed in previous papers (References 2 and 
3). 


In this paper attention will be focused on trajectories that are launched from 
the surface of the Earth and employ a parking orbit before being injected into a 
translunar trajectory. It will also be restricted to an elliptical Earth to Moon 
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TO MOON 



LUNAR TRAJECTORY 


(a) Lunar injection from a direct launch 


TO MOON 



(b) Lunar injection from a parking orbit 

Figure 1 -Trajectories for a lunar mission 
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transfer trajectory, with respect to the Earth, which intersects the Moon before 
the apo -apsis. 

D. Parking Orbit Consideration 

It is envisioned that capability must be provided for the parking orbit to in- 
clude more than one revolution around the Earth. The purpose of this is to pro- 
vide time to check out the various systems on-board the spacecraft and to estab- 
lish the orbit elements so the position, time and velocity necessary for sending 
the spacecraft to the Moon can be calculated. In general, the spacecraft may be 
injected to its translunar trajectory once per revolution from its parking orbit. 


IV. SOLUTION OF THE PROBLEM 

A. Time of Launch into Proper Plane 

As stated previously, one of the requirements that must be satisfied for a 
vehicle leaving the surface of the Earth and arriving in the vicinity of the Moon 
with a minimum amount of energy expended is that the vehicle be launched into 
a plane that is common to the launch site at the time of launch and the Moon at 
the time of impact. This plane may be defined if the declination of the launch 
site, the launch azimuth, the position of the Moon at the desired day and time of 
impact are given. In addition, the times of launch per day into this plane may be 
found if the longitude of the launch site and the desired day of launch are known. 
In this section the equations for finding the time of launch will be presented. 

Consider the geocentric equatorial coordinate system, x , y, and z shown 
in Figure 2. Associated with these coordinates are the unit vectors I , J, andK. 

The unit vector r ; is directed toward the launch site at time of launch, t, , 
on the day of launch, D; , and the unit vector r is directed toward the Moon at 
time of impact, t , on the day of impact, D • The plane that contains both unit 
vectors has an inclination angle, i, at the equator and an azimuth angle, A / , at 
the launch site. The launch site is defined by the right ascension, a l , and the 
declination, S 7 . The unit vector w is normal to the plane defined by t 1 and ? m . 

The following equations may be written 


w z = cos i = cos Sj sin A ; 


( 1 ) 
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Figure 2-Coordinate system 


9 


which is derived in Appendix A 


w • r = 0 


( 2 ) 


w • w = 1 


(3) 


where 


w = I w + Tw + Kw, 

V ^ V Z 


(4) 


and 


r m = If mx + J r n,y + Kr mz 


(5) 


Equations (2) and (3) may be solved simultaneously (Appendix A) to yield 


w = 
y 


-w r r + r 

z my mz — mx 


Vl - r 2 - w 2 

r mz z 


( 6 ) 


r 2 + r 2 

mx my 


— w r - w r 

w y my z mz (7) 


Examining Equations (1), (6), and (7), it is seen that the launch declination, 
launch azimuth and position of the Moon at impact define two planes that are 
fixed in space. Launch into these planes is possible when the launch site on the 
rotating earth passes through the planes and the launch azimuth is in the plane. 
This occurs, in general, twice each day. Figure 3 is included to further clarify 
the geometry of the problem. Launch will be possible once a day or not at all 
depending on whether the value of the radical in Equation (6) is zero or imaginary, 
respectively. 

Referring to Figure 2, additional equations that may be written are: 


• w = 0 


(8) 


and 


(w x r ; ) • K = cos 7 ] = cos 8 ; cos Aj 


(9) 
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LAUNCH SITE AT 
FIRST LAUNCH 



which is derived in Appendix A where 


r ; = I cos a t cos S + J sin cos S ; + K sin 8j 
Equations (8) and (9) may be solved simultaneously (Appendix A) to yield 


cos a 


w w sin 8, + w cos 8, cos A. 

z x l y l l 


l 


(ID 


(w^ - 1) cos 8j 


sin a t 


w w sin 8, - w cos 8, cos A, 

y z /x l l 

(w l - 1) COS S, 


( 12 ) 


and the right ascension, a ; , at the time of launch may be found. Figure 4 shows 
the relationships that exist between the time and angles on the day of launch. 

The following equation may be written for the time of launch on the launch day. 


a { - V - L ; 


(13) 


where is the Earth's absolute rate of spin per mean solar hour, the angle V 
is the mean sidereal time at Greenwich, given by the expression 

t D ' 

¥ = 1.72218633 + 1.720279168 x 10“ 2 D, + .67558729 x 1 0" 5 ( 

where ¥ is in radians and D, is the number of days from January 0.0 U.T. 1960 
to 0.0 U.T. on the day of launch. 

Equation (14) is Newcomb's expression given in Reference (5) with the units 
changed to radian measure and the reference time changed from Noon 1900 Jan- 
uary 0 at the Greenwich meridian to Midnight 1960 January 0 at the Greenwich 
meridian (Appendix A). 

Since the time of lunar impact was assumed and the time of launch calculated, 
the total time of the flight in hours may be obtained 

T « = *. - ‘i + 24 <0„ - D,) . T„ (15) 
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EARTH 



(ZERO HRS. UT) 

a l - + ajt ; +L ; 


Figure 4-Position of the launch site relative to Greenwich 


In summary, the time of launch for a vehicle going from the surface of the 
Earth to the Moon without a plane change may be found if the launch declination, 
longitude and azimuth at the time of launch as well as the position of the Moon 
at the time of impact are known. 

B. Motion in the Plane 


Once the launch time is established, the vehicle can be launched into the 
proper plane that will intersect the Moon at the time of impact. Launch into the 
proper plane is a necessary but not a sufficient requirement for the vehicle to 
impact the Moon at the preselected time of impact. Figure 5 shows the geometry 
in the space fixed plane in which the vehicle is moving after launch. At the pre- 
selected time of impact with the Moon, t m , the Moon passes through the plane at 
the position denoted by the vector R . In order to impact the Moon, the vehicle 
must arrive at the proper position in the plane at exactly the right moment. 

Since a parking orbit approach is to be used in this paper, the problem of the 
in plane motion may be broken into four phases: the first boost, the parking 
orbit, the second boost, and the translunar trajectory. The conditions for a lunar 
impact are then that the sum of the times of each of these phases must be equal 
to the total time given by Equation (15) and the vehicle must be at the position 
denoted by the vector R m . 

1. First Boost — The first boost phase of the problem begins at the time of 
launch and terminates at injection into a circular parking orbit. Both the total 
boost time, T bl , and the boost angle, £ bj , will be assumed known constants. The 
total first boost time may be expressed by the following equation: 


2. Parking Orbit — The parking orbit phase of the problem begins at time of 
injection into the parking orbit (end of first boost) and terminates at the beginning 
of the second boost. Since a circular parking orbit will be assumed, the total 
time in the parking orbit can be expressed by the equation 


Tpo ~ *2 *1 ' ^1 ^po 


(17) 


where C x 
given by 


is the reciprocal of the orbital rate of the parking orbit (Reference 6) 


C ! = 



(18) 
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Figure 5-Geometry of the transfer ellipse 


In Equation (18) /j. is the Earth’s gravitational parameter, Rj is the known 
magnitude of the parking orbit vector, and C po is the parking orbit angle, ex- 
pressed by (See Figure 5) 


£po =( - C b . - - £ f 


(19) 


in which £ b( and £ b2 are known constants, and the angular distance, £, from 
launch to impact may be obtained from the following equation 


f = cos 1 (r, • r m ) 


( 20 ) 


which may be rewritten 


<£ = cos -1 (r cos a, cos S, + r sin a, cos 8. + r sin b,) (^) 

b v mx / l my / l mz / 

The translunar trajectory angle, £ f , is obtained from the expressions given 
in Section (4) following. 

3. Second Boost - The second boost begins at the time the vehicle leaves 
the circular parking orbit, t 2 , and terminates at injection into the translunar 
trajectory, tj . The total boost time, T b2 , and the boost angle, £b 2 » an( * injection 
radius, Ri , will be assumed known constants. The total boost time may be ex- 
pressed by the following equation. 


2 


= t. 
1 


- t„ 


( 22 ) 


4. Translunar Trajectory — The translunar trajectory is the portion of the 
flight between the injection time, t. , and the time of lunar impact, t m . The total 
time of flight during the translunar trajectory, T f , is given by the expression 


In order to study this portion of the flight, it was assumed that the translunar 
trajectory could be approximated by the two-body equations which neglect the 
effect of the Moon gravity on the vehicle. References (2) through (4) indicate that 
this assumption is adequate in obtaining first order estimates of injection condi- 
tions. Utilizing the above assumption, Equations (24) through (31) may be written 
(Reference 6). These lead to a solution of Equation (23). 
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( 24 ) 



where the mean anomaly, M , is given by Kepler's equation 


M. = E. - e sin E. ; j = m, i 
the eccentric anomaly, E, may be obtained from 

E j =cos_ 1 (“ir i ) ; i =m - i 

where 


(25) 


(26) 



(27) 


and 

p = R. (1+ecosv. );j=mori (28) 

The expressions necessary to find the true anomaly, v, and the eccentricity, 
e , are given by 


v m - cos -1 (2RV 2 cos 2 y. - ljJ (29) 


v. = cos -1 (2 V 2 cos 2 y. - l)J (20) 

and 

r ^ -.1/2 

e = I 4 .V 2 (V 2 - 1) cos 2 y. + ll (31) 

The particular form of these equations are derived from their more familiar 
form in Appendix A. 

Upon assuming an injection velocity ratio, V, for an injection elevation angle, 
y i , Equations (29), (30), and (31) are solved which in turn yields a solution for 
Equation (23). 
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The value of the translunar trajectory angle, £ f , necessary to solve Equation 
(19) may be found from 


£ f 


v 

m 


V . 

1 


(32) 


The total time of the flight, T t , can be expressed by the equation 


T = T 

t A b, 


po 


tB 


(33) 


This equation can be solved since T bl , T po , T b2 and Tf are known from Equations 
(16), (17), (22), and (23) respectively. 


It will be noted that the total time of flight, Tf , had been previously calculated 
by Equation (15). Therefore, for a solution of the problem to exist, the total time 
of flight, Tf , in Equations (15) and (33) must be equal. If they are not equal, then 
another velocity ratio, V, is assumed; and the total flight times calculated again. 
This iteration process is continued until the total flight times in Equations (15) 
and (33) are equivalent. At this time the velocity ratio which yields a solution to 
the problem, V f , has been found. 

5. Range of Velocity Ratios and Flight Times — As previously stated, the 
translunar trajectory will be an ellipse, with respect to the Earth, which intersects 
the Moon before the apo-apsis. This type of trajectory yields a specific upper 
and lower limit between which the velocity ratio, V, is to be assumed. 

The lower limit, V x , may be calculated by the following equation (derived 
in Appendix A). 


The upper limit, V 2 


V i = 


1 - R 


-I 1/2 


I_l - (R cos Ti) 2 J 


is taken to be that of parabolic flight. 


(34) 


V 2 - 1.0 (35) 

The velocity ratio which yields a solution to the problem, V f , lies between 
these two limits. In order to get an idea of the values of these limits, Equation 
(34) may be solved for its minimum anticipated value. To do this some numbers 
consistent with the ideas previously proposed must be assumed. 
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Assume: 


a. The earth to moon distance, , varies between a maximum of 63.8 
earth radii to a minimum of 55.8 earth radii. These values were ob- 
tained from Reference 7. 

b. The earth injection distance, Rf , varies between 1.0 earth radii and 1.1 
earth radii. These values seem adequate to cover the injection range of 
current boosters. 

c. The injection elevation angle, 7i , varies between zero degrees and 20 
degrees. This seems to be a reasonable assumption. 

Using these assumptions, Equation (34) has been evaluated and the results tabu- 
lated in Table I. From this table it is seen that the minimum anticipated velocity 
ratio, min V x , is 0.990264. The range of velocity ratios that will yield a solution 
is, therefore, seen to lie between 0.990264 and 1.0. 

Since one of the given pieces of data will be the day of launch, D ; , it is im- 
portant to find what the range of flight times from injection to lunar impact, T f , 
are. This may be done by utilizing the previous assumptions and corresponding 
velocity ratios. To calculate the maximum anticipated flight time, max T f , 
Equations (23) to (31) must be evaluated using V x . 

This has been done and the results appear in Table I. From this table it is 
seen that the maximum anticipated flight time, max T, , is 130.1477 hours. To 
calculate the minimum anticipated flight time, min Tf , the value of V 2 must be 
used. Since this is a parabolic flight, Equations (24) to (28) are no longer valid. 
Instead they are replaced by the following equations (Reference 6). 


T , = + R » + s > 3/! - < R i + R . - s > 3/3 ] 


(36) 


which holds for parabolic flight when v m - v. < -n and where 


S = [r 2 +R 2 - 2R.R cos (v - i/.)l 

[_ i m im v m i J 


1/2 


(37) 


The calculation is then made using Equation (35), Equations (29) to (31), and 
Equations (36) and (37). Results appear in Table I where it is seen that the 
minimum anticipated flight time, min T f , is 44.9697 hours. The range of flight 


19 



Table I 


Range of Velocity Ratios and Flight Times 


Assumed 

Calculated 

R 

m 

R 

1 

7i 

v. 


v x 

T f 

earth radii 

earth radii 

degrees 

dimensionless 

dimensionless 

hours 




0 

* 


.991158 

106.5595 

55 

.8 

1.0 

20 

4 


.991139 

106.1230 




0 

* 


.990287 

106.8411 



1.1 

20 

* 


.990264 

106.3542 

1 

1 


0 

* 


.992253 

129.8469 

63 

.8 

1.0 

20 

* 

.992240 

129.3886 




0 

* 

.991489 

130.1477 



1.1 

20 

* 


.991472 

129.6370 

'I 

r 


0 

1.0 

* 

45.2046 

55 

.8 

1.0 

20 



* 

44.9697 




0 



* 

45.3196 



1.1 

20 



* 

45.0562 


f 


0 



* 

55.0896 

63 

.8 

1.0 

20 



* 

54.8449 




0 



* 

55.2131 

1 

r 

1.1 

20 

1 


* 

54.9389 


times that will yield solutions is, therefore, seen to lie between 130.1477 and 
44.9697 hours. Knowing this, along with the desired day of lunar impact D ra , and 
the desired number of parking orbit revolutions, it is seen that solutions of the 
problem will exist anywhere between seven and one days before lunar impact. 
Therefore, a good choice for D, to obtain all solutions would be in the range from 
D - 7 to D - 1. 

m m 

6. Obtaining Successive Velocity Ratios — Having established a lower limit 
velocity ratio, , and an upper limit velocity ratio, V 2 there remains a pro- 
cedure to be followed in obtaining successive assumed velocity ratios that will 
ultimately lead to a solution. There are many methods by which this iteration 
may be done, the one described below is known as the Regula Falsi method of 
iteration. Calling T t from Equation (15) T tA and T t from Equation (33) T tB , 
the quantity AT t may be calculated by 


(AT t ) k 


<?t B \ - T tA 


k = 1, 2, 


(38) 


20 








The object of the iteration is to yield a AT t , which is a function of V, equal to 
zero. Figure 6 shows the procedure when AT t is plotted as a function of V. 

After Equation (34) is used in calculating Vj , this value is used along with other 
knowns in calculating T tB , called (T tB )j since Vj was used in its calculation, by 
use of Equations (23) to (33). Since T tA is already assumed known, and constant 
for this procedure; the maximum value of AT t) viz. (AT t )i , may be calculated 
by Equation (38). The upper limit, V 2 , is known and another value of T tB , viz. 
(TtB) 2 > is calculated by use of Equations (29) to (33) along with Equations (36) 
and (37), since the trajectory is parabolic. This quantity along with the constant 
and known T tA is used in evaluating Equation (38) for the minimum value of AT t , 
viz. (AT t ) 2 . Through the two points on the solution curve, [(ATt)!, Vj] and 
[(AT t ) 2 , V 2 ] , a straight line (straight line 1) is constructed. The next value of 
the velocity ratio to be assumed, V 3 , is located at the intersection of straight 
line 1 and AT t = 0. Analytically the expression yielding this new assumed velocity 
ratio and successive assumed velocity ratios is found to be 


'k + l 


(AT,), V„ - (AT,) k V, 
(AT ) - (AT,) 

1 1 k 


; k = 2, 3, 


(39) 


The value of V 3 is then used to calculate (T tB ) 3 by use of Equations (23) to (33). 
A new minimum value of AT t , viz. (AT t ) 3 , may then be calculated by Equation 
(38). A new straight line (straight line 2) is constructed through the two points, 
[(^Tt),, Vj] and [(AT t ) 3 , vj, on the solution curve. The next velocity ratio to 
be assumed, V 4 , is the intersection of straight line 2 and AT t = 0. Analytically 
this also may be found by Equation (39). The iteration continues in this manner 
yielding values of AT t closer and closer to zero as it proceeds. When (AT t ) k 
equals zero, the point where the solution curve intersects AT t = 0 has been 
found, and the velocity ratio yielding a solution to the problem, V f , is known. 

7. Initial Injection Conditions — If a solution to the problem exists, a value 
for v f is found; and it is then possible to calculate the initial injection conditions 
for a translunar trajectory. These conditions can be obtained in either cartesian 
coordinates (x. , y. , z , x. , y , z. ) or polar coordinates (L. , S. , R. , V. , A. , 

The cartesian coordinates may be found by the following equations. 


X. 


R. 

1 


1 


= E, ( -a,) R 2 ( 8, ) R, (A, - n/2) R, ( C po -C b; ) 

0 

z . 

i 


0 


% 
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SOLUTION CURVE 
— STRAIGHT LINE 1 


STRAIGHT LINE 2 




and 


z. 

1 


R 3 (-a,) W V A /--/ 2) R 3 (-S bi -^-^-^2+^) 


where the rotations R x (0), R 2 (^), andR 3 ( 9 ) are defined as 


V*) 


10 0 
0 cos 9 s in 6 


R,(0) = 


R, ( 0 ) 


(41) 


(42) 


[_0 -sin# 

cos 8 


cos 0 

0 

-sin 8 


0 

1 

0 

(43) 

sin# 

0 

cos 8 


cos 6 

sin 9 0 


-sin 8 

cos 9 0 

(44) 

0 


0 1 _ 



^ The injection velocity, V. , may be obtained since the final velocity ratio, 

V f , is known, and the parabolic velocity at injection, V p , is given (Reference 6) as 


Therefore 



(45) 


V. 



(46) 
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In polar coordinates three of the quantities Ri , V, , and Ji are known by 
products of the solution, and the remaining quantities may be found by the equations 


L. 

1 



- (y + a)t l ) - cu(t { 


*/> 


(47) 


and 



(48) 


A. = tan 1 
1 


R. (x.y. - y.x. ) 

1 V l* 7 1 J 1 1 ' 


[y. (y.i. - z.y.) -x. (z.x. - x.^jj 


(49) 


8. Solutions for Successive Orbits - The problem as outlined thus far will 
yield a solution for the first revolution of the parking orbit around the earth. To 
obtain a solution for successive orbits, the following technique is used. 

Figure 7 is a sketch of the geometry for the solutions in the plane of the 
trajectory. The position of injection for the first parking orbit solution is indi- 
cated by the vector (R A ) j and the position of injection for the second parking 
orbit solution by vector (R ) 2 . Note that the second position occurs after re- 
maining in the parking orbit somewhat more than one full revolution in inertial 
space. This comes about from a combination of facts. 

a. The time of launch from the surface of the earth to the time of lunar 
impact remains the same. 

b. The injection elevation angle, y i , is the same. 

c. The moon is essentially in the same position since the time for one 
parking orbit 1.5 hrs.) does not allow for the moon to move very 
much (moon's motion 12°/day). 

d. Since more time is spent in the earth parking orbit as the number of 
parking orbits is increased, it is evident that the time in the translunar 
trajectory must decrease in order for a. to be true. 

The above facts necessitate a higher energy translunar trajectory which 
goes through essentially the same point indicated by the vector R m . The only 


24 



PARKING ORBIT 
ALTITUDE 



way this can be accomplished is by rotating the position injection vector from 
(R i ) x to (R i ) 2 in a counter-clockwise sense as shown in Figure 7. 

To find solutions for successive orbits, the procedure is to subtract (n - 1) 
2 77 from the total time given by Equation (15) or 


<V„ = *.-«,+ 24 (D. - D,) - (n-1) 2 „ C, - <T tJ )„ < 5( 

where n is the number of the orbit in which the solution is to be found. In all 
other respects the problem remains the same. 
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Appendix A 
DERIVATIONS 


1. Derivation of Equation (1): Refer to Figure Al, and use the law of sines 


s in (90° - i) _ sin (90° - Sf) 
s i n A ; sin 90° 


w = cos i = cos S, sin A, (^) 

2 ll 

2. Simultaneous solution of Equations (2) and (3) to yield Equations (6) and (7): 

w • r = 0 ( 2 ) 

m 

w • w = 1 (3) 


w • r = w r + w r +wr =0 

m x mx y my z mz 


W ‘ W = W 2 + W 2 + W 2 = 1 


w = 

x 


-(w r + w r ) 

y my z mz ' 


( 7 ) 


W 2 — 1 — W 2 — W 2 — 1 — 
y x z 


' w r + w r \ 46 

y my z mz \ 7 

r - W ' 

mx / 


Which upon simplification becomes 


(r 2 , r 2 ) w 2 + (2w r r ) w + (w 2 r 2 +w 2 r 2 -r 2 ) = 0 

' mx + my' y ^ ' z my mz ' y v z mx z mz mx' 
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Solving this quadratic equation in w y yields 


-w r r ± r Vi - r 2 - w 2 

z my mz mx ¥ mz z 

W = 

y 2 2 
r 1 + r ^ 

mx my 

3. Derivation of Equation (9): Refer to Figure A2, it is seen that 

(w X Fj ) • K = COS Tj 

The direction cosines of K yield 

cos 2 i + cos 2 77 + cos 2 (90 - Sj ) = (K) 2 = 1 
Substituting in Equation (1) for the first term, transposing, and simplifying 

cos 2 rj = 1 - cos 2 Sj sin 2 Aj - sin 2 8 j = cos 2 Sj cos 2 Aj 
Upon taking the square root of both sides 


COS Tj - cos S i cos Aj 

where the positive value is considered only if launch is in the northern hemisphere 
(n/2 > Sj > 0) and eastward ( tt > > 0) 


/. (w x r^) • K = cos 77 = cos 8; cos Aj 


(9) 


4. Simultaneous solution of Equations (8) and (9) yield Equations (11) and (12). 


? i * w = 0 


( 8 ) 


(w x r t ) • K = cos cos A j 


(9) 


Since 


r t - I cos ctj cos 8, + J sin a j cos 8; + K sin 8 ; 


( 10 ) 
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Performing the vector operations yields 


r l • w = cos aj cos S ; w x + sinaj cos Sj w + sin w z = 0 


(Al) 


(w x ) • K = sin a t cos S l w x - cos a, [ cos 8; w = cos 8 [ cos A ; (A2) 


Upon substituting (A2) into the second term of (Al) 

w 

cos a l cos w x + — (cos a l cos Sj w y + cos S ; cos A ; ) + sin S ; w z = 0 

W x 

Which upon transposing and solving for cos a l yields 


w w sin 8, + w cos 5. cos A, 

z x i y * l 


cos a 


(ii) 


(w^ - 1) cos 8 j 

Upon substituting (Al) into the second term of (A2) 


w 

sin cos 5^ w x - (-sin cos Sj w y - sin Sj w z ) ~ cos Sj cos Aj 

w x 

which upon transposing and solving for sin a l yields 


sm a 


w w sin 8, - w cos 8, cos A » 

y Z lx l l 


( 12 ) 


(w 2 - 1) cos 8j 


5. Ch ang ing the reference time and units of measure in Newcomb's expression. 

Greenwich mean siderial time at zero hours universal time on successive 
dates are computed from Newcomb's expression (Reference 5). 


¥ = 6 h 38 m 45?836 + 8640184?542 T + 0?0929 T 2 (A3) 

O 

where for any date, T denotes the number of Julian centuries of 36525 days which 
have elapsed since noon on 1900 January 0.0 at the Greenwich meridian. 
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This equation may be written as 


1^hoo+ AT + b T : 


To shift the reference time from January 0.5, 1900 to January 0.0, 1960 


*P 0 = ^1900 + A ( T + T i) + B ( T + V 2 = figoo + AT + bt2 + (A + 2 BT) Tj + BT 2 


where T is the time in Julian centuries between January 0.5, 1900 and January 
0.0, 1960 and Tj is the time in Julian centuries after January 0.0, 1960. The 
value of T in this case is 


T = x ^0 + 13-5 _ 599950932 Julian centuries 
36525 


*o = ^1 9 6 0 + t (8640184.542 + 2 ( .0929) ( .599958932)] T x + .0929 T 2 
and the Equation (A3) becomes 

W o = 6 h 34 m 41? 762 + 86401 84? 653 T x + 0?0929 T 2 (A4) 

where T x denotes the number of Julian centuries which have elapsed since 1960 
January 0.0. 

The units of Equation (A4) may be changed to yield 


¥ = 1.72218633 + 1.720279168 x 10" 2 D, + 0.67558729 x 10 “ 5 f — ^ — ) ( 14 ) 

* \36525 / 

where V is in radians and Dj is the number of days from January 0.0 U.T. 1960. 

6. Equations (29), (30), and (31) may be obtained in the following manner. Start- 
ing with the basic equations (Reference 6). 

p = cos 2 y (A5) 
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p = a(l-e 2 ) 



Equating p in Equations (A5) and (A6) and solve for e 2 



Substitute for l/a using Equation (A7) 

** - - [ 4 -?] [■©'“•'>] 

Equation (A7) also implies for parabolic velocity (i.e. a =co). 



Substitute this relationship in and solve for e yields 

e = ^4 V 2 (v 2 - l) cos 2 7. + lj 


(A6) 


(A7) 


(A8) 


(A9) 


(31) 
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Combining Equations (A5), (A9), and (A8) yields 


COS V 


, = r (2 v 


2 cos 2 y - l) 


(30) 


Applying Equation (A5) for earth injection distances and lunar distances yields 
the results 


P 



cos 2 y 

/ m 



cos 2 y. 

' i 


Substituting? /R m from above results along with (A9) into Equation (A8) yields 


cos 


- 1 


- (2 R V 


2 cos 2 y. 


- 1 ) 


] 


(29) 


7. Derivation of the lower limit for assumed velocity ratios. 

Equation (35) for the lower limit, V , is obtained in the following manner. 
Upon substitution of Equation (31) into Equation (29) and squaring yields 


cos 2 v 

m 


4 R 2 V 4 cos 2 y. - 4 R V 2 + 1 
4 v 2 (v 2 - l) c o s 2 y. + 1 


(A10) 


The limiting case occurs when cos 2 v m = 1 which corresponds to a lunar 
impact at the apo-apsis of the translunar ellipse. Making this substitution and 
solving for the velocity ratio Equation (A10) becomes 


1 - R 


-|l/2 


(_ 1 - (R cos y.) 2 J 


(34) 
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Figure A-l -Geometry of the orbit 
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Figure A-2-Geometry describing the angle 7 ] 
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APPENDIX B 


CALCULATING PROCEDURE 


A. Data Required 

In order to solve the problem of a lunar impact from a given launch site, 
various data are required before the calculation may begin. A summary of these 
input data and its source are given below. 

1. The first boost are £ b i : Refer to Figure 5. This is a characteristic of 
the booster to be used in the mission. It is usually available or may be 
obtained from data in reports giving the booster characteristics. 

2. The first boost time, T bl : Refer to Figure 5. The same remarks in 
Item 1 apply here. 

3. Parking orbit altitude, h po = R t - R e . The altitude at which the first 
boost injects its payload into a circular orbit. The same remarks in 
Item 1 apply here. 

4. The second boost arc, £b 2 : Refer to Figure 5. This is a characteristic 
of the booster used to boost the vehicle from circular parking velocity 

to translunar injection velocities. This information may be obtained from 
data in reports giving the booster characteristics. 

5. The second boost time, T b2 : Refer to Figure 5. The same remarks in 
Item 4 apply here. 

6. Injection altitude, h* = R t - R e . The altitude at which the second boost 
injects the vehicle into its translunar trajectory. The same remarks in 
Item 4 apply here. 

7. Injection elevation angle, y. : Refer to Figure 5. This angle is a function 
of the particular guidance system employed during the second boost. 

8. Launch azimuth, . The desired launch azimuth at launch site. Depends 
upon range safety and tracking requirements. 

9. Parking orbit revolutions. The number of parking orbit revolutions for 
which solutions are desired. 
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10. Launch days(s), D 7 . The number of days before impact that the launch 
is to take place. For the elliptical trajectories described this is from 
7 to 1 days before D m . 

Flight time, T f , the time of the flight from injection to impact may be 
given instead of D, . 

11. Lunar impact time, Y ra , D m , t m . The year, day and time it is desired 
to impact the moon. 


B. Calculating Steps 

A flow diagram indicating the solution of the problem is given in Figure Bl. 
This figure is a guide to the solution of the problem and indicates the general 
flow from equation to equation in the text. Steps in the calculating procedure 
which may be used in conjunction with this figure are given as follows: 

1. Given a year, Y ra , day, D m , and time, t m , it is desired to impact the Moon. 
The position of the Moon in geocentric equatorial coordinates, R^, Rmy , 
R mz may be obtained from an ephemeris of the Moon. The resultant of 
these coordinates, !?„ , as well as the components of the unit vector, r mx , 
r my ’ r mz > are then calculated. 

2. Given the launch site declination, S / , and a launch azimuth, A / , compon- 
ents, w x , w y , w z of the unit vector w normal to and defining the plane of 
motion from the surface of the Earth to the Moon are obtained from Equa- 
tions (1), (6), and (7). Depending upon the value under the radical in 
Equation (6) the unit vector w will have two values, one value, or be imag- 
inary. These results indicate for the launch azimuth, A; , launch into the 
plane is possible twice a day, once a day, or not at all. If a real value 

of w exists the calculating procedure is continued. If two real values of 
w exist, use one of them at a time. 

3. The right ascension at launch, a ; , may be calculated by Equations (11) 
and (12). A desired day of launch, D/ , is given (i.e. D / is an integer in 
the range of D m - 7 to D m - 1) and the angle is calculated by Equation 
(14). Since the launch longitude, L; , and the Earth's spin rate, oj , are 
knowns, the time of launch after zero hour U.T. on the day of launch, D; , 
can be calculated from Equation (13). The day of launch, D/ , and the time 
of launch, t f , as well as the day of impact, D ra , and time of impact, t m , 
are now all knowns; the total time of the flight, T t , can be calculated from 
Equation (15). This time is designated as T tA . 
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4. A value of the injection angle, 7i , is given, and the injection radius, Ri , 
can be calculated from the given injection altitude, hi , and Earth's radius 
R e . A velocity ratio, Vj , calculated by Equation (34) is now obtained. 

5. Using this velocity ratio, the eccentricity, e , as well as the true anomalies 
to the Moon, v m , and to the injection point, , are calculated by Equations 
(31), (29), and (30) respectively. The translunar trajectory angle, , 
may now be calculated from Equation (32); and the angular distance, f , is 
calculated by Equation (21). Since the booster burning angles, £b, and 

£b 2 > are given, the parking orbit angle, £ po , can be calculated from Equa- 
tion (19). 

A circular parking orbit was considered in this analysis; and therefore, 
the parking orbit radius, Ri , is constant and is obtained from the given 
parking orbit altitude, h po , and Earth's radius. The inverse orbital rate 
of the parking orbit, Cj , is defined by Equation (18). This is used in 
calculating the time in parking orbit, T po , given by Equation (17). 

6. Since the velocity ratio yields an elliptical flight path, the calculation 
continues as follows: The semilatus rectum, p , and the semimajor 
axis, a , may be calculated by Equations (28) and (27) respectively. Hav- 
ing these quantities the eccentric anomalies to the Moon, Em , and to the 
injection point, Ej , can be calculated by Equation (26); and the mean 
anomalies to the Moon, M m , and to the injection point, Mi , can be calculated 
by Kepler's equation, viz. Equation (25). The time difference on the trans- 
lunar trajectory between injection and lunar impact, t m - tf t can be cal- 
culated by Equation (24). This time difference, called T f , is the same as 
that defined by Equation (23). 

7. The boost times, Tb, and T b2 , are both given and are used with T po and 
T f to calculate the total time of the flight, T t , from Equation (33). This 
time is designated as T tB . In step (3) above, it is seen that the total time 
of flight, T t = T tA , has already been calculated by another method, viz. 
Equation (15). The difference in these times, AT t , is given by Equation 
(38). 


8. The particular AT t as calculated by Equation (38) for Vf from Equation 
(34) is designated, since k = 1, as (AT t )j on Figure 6. If (AT t )j is less 


36 



than zero, the solution of the problem requires that the transfer trajectory 
intersect the Moon after apo-apsis and the calculation is concluded. A 
new value of the launch day, D; , greater than before is used, and the cal- 
culation is reinitiated at step (3), Equation (14). If (AT t ) j is greater than 
zero, continue. 


9. The next velocity ratio used, V 2 , is given by Equation (35). Step (5) is 
then repeated. Since the velocity ratio corresponds to a parabolic trajec- 
tory, the parameter, S , is calculated by Equation (37) and is used in cal- 
culating the total time from injection to lunar impact, T f , given by Equa- 
tion (36). Step (7) is then repeated. The particular AT t as calculated by 
Equation (38) for V 2 from Equation (35) is designated, since k = 2, as 
(AT t ) 2 on Figure 6. If (AT t ) 2 is greater than zero, the solution of the 
problem requires a velocity ratio greater than parabolic; and the calcu- 
lation is concluded. A new value of the launch day, D/ , less than before 
is used, and the calculation is reinitiated at step (3), Equation (14). If 
(AT t ) 2 is less than zero, a new velocity ratio, V 3 , must be calculated. 

10. The next velocity ratio to be assumed may be calculated by Equation (39). 
Steps (5), (6), (7), and (10) are repeated in that order and are continued 
to be repeated until the value of AT t given by Equation (38) is as close to 
zero as desired. When this occurs, a velocity ratio which has yielded a 
solution to the problem, Vf , has been obtained. 

11. At this point several things may be calculated depending upon the partic- 
ular purpose of making the calculation. Two possibilities come to mind, 
and to distinguish them one will be called Option A and the other Option B. 

Option A: In this case it is assumed the purpose for making the cal- 
culation is to obtain all the possible injection conditions that will yield a 
lunar impact at a given desired time for given launch azimuth, booster 
conditions, and injection elevation angle if the parking orbit is allowed to 
orbit the earth several (sayn ) times. For this option n was assumed to 
be equal to one through steps 10. The value of n is now increased to two, 
and Equation (50) is used in place of Equation (15), the entire calculation 
being repeated from step 4. This value of n is increased as many times 
as desired repeating the calculation from step 4. Having done this, a 
solution has been found for each revolution of the parking orbit (i.e. n 
of them). Assuming two values of w y were obtained when Equation (6) 
was solved (see step 2), there remains now to repeat the calculation from 
steps (3) to (10) for the second possible value of w y . Having done this 
(for n = 1), the values of n are increased again one at a time as was 
previously done. Now the two possible solutions per day for each parking 
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orbit revolution have been found, yielding 2n possible solutions. The 
results thus far have been for only one assumed launch day, D ; . Since 
lunar impact at the desired time is possible, as previously explained, 
for launch times between 44.9697 and 130.1477 hours before lunar im- 
pact, more launch days must be investigated. To be sure to include all 
possible cases, the launch days to be investigated are between D m - 7 
and D m - 1, or a total of 6 days. Doing this for n revolutions in the 
parking orbit for each day is seen to yield a possible 12 n solutions per 
launch azimuth per elevation angle that will impact the moon at the same 
time. 

Option B: In this case it is assumed the purpose for making the cal- 
culation is to obtain initial injection conditions that will have a given 
desired flight time, T f , from injection to lunar impact. Lunar impact 
will be near the desired time, but not at it. Solutions will be obtained 
in the first parking orbit for a given launch azimuth, booster conditions, 
and injection elevation angle. The calculation is made as indicated in 
steps 1 through 10, and a solution is obtained. At this point the desired 
flight time (Tf ) D , is compared with flight time, Tf , that was the solution 
to the problem, and their difference, AT f = (ATf) D - T f , is obtained. 

This time difference is algebraically added to the original lunar impact 
time (thus having the effect of moving the moon), and steps 1 through 10 
are again repeated, yielding a new solution. Again a time difference 
AT f , is obtained, and the process is repeated until this time difference 
is zero. A solution during the first revolution of the parking orbit having 
a given time of flight from injection to impact has thus been found. This 
solution will yield a lunar impact near (say within a day) the original 
impact time desired. The original impact time is then utilized with the 
second value of w y obtained when Equation (6) was solved, and the entire 
process of Option B is repeated. The final results will yield a total of 
two solutions, having the same desired flight time, in the first parking 
orbit for a given launch azimuth, injection elevation angle, and boost 
conditions that impact the moon close to the desired impact time. 

12. Having found the final velocity ratio, V f , all the by-products necessary 
for finding the injection vector components in both cartesian coordinates 
( x i. y i > z i, x i , Yi , z \ ) or polar coordinates (Lf, S, , Ri , V, , ,7;) 

may be found. Parabolic velocity, V p , is found using Equation (45) which 
then enables the calculation of Vj from Vf, Equation (46). All the neces- 
sary data is now known to calculate rotation matrices, Equations (42), 
(43), and (44), which are used in obtaining the cartesian position com- 
ponents and velocity components from Equations (40) and (41) respec- 
tively. Using these components the remaining unknown polar coordinate 
components, viz.Lj , , and Aj , may be calculated by use of Equation 

(47), (48), and (49) respectively. 
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Note: Numbers ia parenthesis refer 
to equations in the text. 


INPUTS (I) 
































CONSTANTS (C) 

§1 , Lj , w, fi, R e 



Figure B-l-Flow diagram indicating the solution of the problem. 
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